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Q ABSTRACT 

m 

y—{ Context. Globular clusters (GCs) are established emitters of high-energy (HE, 100MeV<E<100GeV) y-ray radiation which could originate 
• • from the cumulative emission of the numerous millisecond pulsars (msPSRs) in the clusters' cores or from inverse Compton (IC) scattering of 
^ relativistic leptons accelerated in the GC environment. These stellar clusters could also constitute a new class of sources in the very-high-energy 
k^l (VHE, E>100 GeV) y-ray regime, judging from the recent detection of a signal from the direction of Terzan 5 with the H.E.S.S. telescope array. 
rS Aims. To search for VHE y-ray sources associated with other GCs, and to put constraints on leptonic emission models, we systematically analyzed 
^ the observations towards 15 GCs taken with the H.E.S.S. array of imaging atmospheric Cherenkov telescopes. 

Methods. We searched for point-like and extended VHE y-ray emission from each GC in our sample and also performed a stacking analysis 
combining the data from all GCs to investigate the hypothesis of a population of faint emitters. Assuming IC emission as the origin of the VHE 
y-ray signal from the direction of Terzan 5, we calculated the expected y-ray flux from each of the 15 GCs, based on their number of millisecond 
pulsars, their optical brightness and the energy density of background photon fields. 

Results. We did not detect significant VHE y-ray emission from any of the 15 GCs in either of the two analyses. Given the uncertainties related to 
the parameter determinations, the obtained flux upper limits allow to rule out the simple IC/msPSR scaling model for NGC 6388 and NGC 7078. 
The upper limits derived from the stacking analyses are factors between 2 and 50 below the flux predicted by the simple leptonic scaling model, 
depending on the assumed source extent and the dominant target photon fields. Therefore, Terzan 5 still remains exceptional among all GCs, as 
the VHE y-ray emission either arises from extra-ordinarily efficient leptonic processes, or from a recent catastrophic event, or is even unrelated to 
the GC itself. 

Key words, globular clusters: general - radiation mechanisms: non-thermal - pulsars: general - gamma rays: general 
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1. Introduction 

Globular clusters (GCs) are very dense systems of stars, hosting 
the oldest and most evolved stellar populations of the Galaxy. 
The extremely high rate of close stellar encounters in the cores 
of GCs likely gives rise to large numbers of dynamically formed 



compact binary systems (Pooley & Hut 2006| |Pooley et al 2003 
Hui et al. 2Q10| > which are believed to be the progenitors of 



millisecond pulsars (msPSRs). Indeed, large msPSR populations 
have been detected in many Galactic GCs through radio obser- 
vations (see Camilo & Rasio 2005; Ransom 2008 ). 

A complementary approach to study the msPSR populations 
of GCs comes from high-energy (HE, 100MeV< E <100GeV) 
y-ray observations. Because msPSRs can be strong sources of 
pulsed magnetospheric emission in this energy regime (as shown 
by observations with Fermi/LAT, Abdo et al.|2009a), y-ray sig- 



nals are expected from GCs either due to bright individual msP- 
SRs or due to the cumulative emission from the whole popu- 
lation. Up to now, 14 Galactic GCs have been detected with 
Fermi/LAT ( fAbdo et al.||2009bl |2010b| |Kong et al.||20T0l [Taml 
et al. 2011) , and the y-ray spectra of most of these sources 
show steep cut-offs at a few GeV, very reminiscent of the spectra 
measured from individual msPSRs. Indeed, an individual y-ray 
msPSR has recently been de tected in the GC NGC6624 with 
Fermi/LAT ( jFreire et al.|20T"fl >. 

Through the observations of very-high-energy (VHE, 
E > 100 GeV) y-rays even the whole msPSR population might 
be accessible. Models predicting emission in the VHE regime 
rely on inverse Compton (IC) up-scattering of soft photon fields 
by ultra-relativistic leptons accelerated in the pulsars' magne- 
tospheres (Venter et al. 2009) or even re-accelerated in col- 
liding pulsar wind nebula shocks (Bednarek & Sitarek 2007). 
Following Cheng et al. (|2010), under certain conditions even the 
HE emission observed by Fermi/LAT could be explained within 
a pure IC scenario. In all these models, depending on the as- 
sumed diffusion timescale, the magnetic field strength and the 
location of the GC, the relevant target photon fields could either 
be the GC starlight, the cosmic microwave background or the 
Galactic optical and infrared background light. In other words, 
if the leptons remain localized close to their sources, e.g. due 
to slow particle diffusion, the photon field of the GC member 
stars would be the dominant target for the IC process. In con- 
trast, if the leptons propagate out of the core region of the GCs, 
other photon fields would become more important, because the 
energy density of the GC photon field rapidly decreases with the 
distance from the core. Thus, in this case the morphology and 
size of the VHE y-ray source would be decoupled from the dis- 
tribution of msPSRs as well as from the extent of the GC stellar 
population. 

Recently, VHE y-ray emission has been detected from the 
direction of the GC Terzan5 with the H.E.S.S. telescope array 
(HESS J1747-248, |Abramowski et al.|201 la) . Should this sig- 



nal indeed be associated to Terzan5 this would establish GCs 
as a new class of sources in VHE y-rays, since previous stud- 
ies of GCs by H.E.S.S., MAGIC, VERITAS, and CANGAROO 
in this energy domain resulted only in upper limits (see e.g. 
Aharonian et al.||2009| |Abramowski et aL"1|2011b[ [A nderhub 
et al.|2009UMcCutcheon for the V ERITAS Collaboration |20091 
Kabuki et al.|2007) . Even though the nature of this source is not 
settled yet, IC emission from leptons accelerated by the particu- 
larly large msPSR population of Terzan 5 remains a possibility. 



In this case also diffuse X-ray emission could be expected aris- 
ing from synchrotron radiation (SR) from the same population 
of leptons (Venter & de Jager 2008). Interestingly, extended and 
diffuse non-thermal X-ray emission was detected from Terzan 5 
with Chandra by Eger et al. ( |2010 i. However, the apparent off- 
set between the X-ray and VHE y-ray signals challenges simple 
(IC/SR) models. A search for diffuse X-ray emission on similar 
scales from other HE-detected GCs resulted only in upper limits 
(Eger & Domainko 2 012) . 

As an alternative to the leptonic msPSR scenarios, a hadronic 
model was proposed for HESS J 1747-248 by |Domainkol ( pOTT] ). 
Here, it is suggested that the VHE y-ray emission arises from 
the decay of neutral pions that were produced in inelastic colli- 
sions of relativistic protons with interstellar matter. These pro- 
tons could have been accelerated in a remnant of a short y-ray 
burst resulting from a merger of binary neutron stars. In the cores 
of GCs the probability for such a rare event is expected to be en- 
hanced compared to Galactic average, due to their large binary 
fraction. W ith a rate for such rare catastrophic events of only 
~10 4 yr _1 |Nakar|2007 1, it is unlikely that two of them can be 
observed at the same time and, thus, Terzan 5 may stand out as 
an exception among all GCs. 

At this point it is thus unclear if or to what extent the X-ray, 
HE and VHE y-ray signals detected from GCs are related and 
what the underlying radiation mechanisms are. Further studies, 
particularly in the X-ray and VHE y-ray regimes are the most 
promising to shed new light on the high-energy processes in 
these systems. In this work we present a systematic search for 
VHE y-ray emission from all Galactic GCs that are covered by 
the present observational data set of the H.E.S.S. telescope ar- 
ray. To increase the sensitivity for the detection of a population 
of individually faint emitters we also performed a stacking anal- 
ysis, where we combined the number of photon candidate counts 
from the individual GCs to test for a statistically significant sig- 
nal. We then compare the obtained results to the VHE y-ray sig- 
nal detected from Terzan 5 to investigate scenarios based on IC 
emission from leptons accelerated by msPSRs. 



2. H.E.S.S. observations 

H.E.S.S. is an array of four imaging atmospheric Cherenkov 
telescopes located in the Khomas Highland of Namibia. Each 
telescope is equipped with a tessellated mirror with a surface 
area of 107 m 2 and a camera comprising 960 photomultiplier 
tubes with a total field of view (FoV) of 5°. For the recon- 
struction of atmospheric air showers the array works in coinci- 
dence mode requiring at least two telescopes to detect individual 
events. In this mode the array features a point-spread-function 
(PSF) with a 68% containment radius of ~0°07 and an energy 
resolution of 15% on average. The stereoscopic reconstruction 
provides an excellent suppression of the cosmic -ray background 
reducing the numbe r of background events after all cuts by a 
factor of up to 10 4 (IBenbow 2005). With modern reconstruc- 
tion techniques the H.E.S.S. array has a point-source sensitivity 
of about 2xl0~ 13 ph cirT 2 s" 1 (E>200 GeV) for 25 hours of ob- 
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servation, requiring a statistical significance of 5 <x (see e.g. de 
|Naurois & Rolland|2009| |Ohm et al.|2009l ). Here a zenith angle 
of 20° and an offset of of 0°7 is assumed. 



2. 1 . The globular cluster sample 



We based our target selection on the GC catalog of |Harris| (2010 
henceforth H10) from which we also took all basic parameters 
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Table 1. The GC sample studied in this work 



GC 


long. (1) 


lat.<» 


zenith (2) 


offset* 2 ' 


livetime' 3 ' 


name 


n 


O 


O 


O 


(h) 


NGC 104 (47Tuc)'*' 


305.89 


-44.89 


49.8 


0.9 


23.1 


NGC6388 W 


345.56 


-6.74 


23.6 


0.7 


17.9 


NGC 7078 (M 15) 


65.01 


-27.31 


37.5 


0.7 


12.3 


Terzan6 (HP 5) 


358.57 


-2.16 


24.6 


1.8 


15.2 


Terzan 10 


4.49 


-1.99 


18.1 


1.6 


4.2 


NGC 6715 (M54) 


5.61 


-14.09 


18.3 


0.7 


11.8 




jOI.jj 


— 40.25 


4y.5 


1 .3 


£ r\ 

3.U 


Pal 6 


2.10 


1.78 


19.4 


1.4 


24.7 


NGC 6256 


347.79 


3.31 


20.5 


1.3 


5.3 


Djorg 2 


2.77 


-2.50 


15.9 


1.7 


4.6 


NGC 6749 


36.20 


-2.21 


33.7 


1.5 


8.2 


NGC 6144 


351.93 


15.70 


26.8 


1.4 


4.7 


NGC 288 


152.30 


-89.38 


11.8 


1.4 


46.7 


HP 1 (BH 229) 


357.44 


2.12 


14.3 


1.5 


5.6 


Terzan 9 


3.61 


-1.99 


18.4 


1.5 


5.2 



'''Galactic coordinates used for the analyses, as given by H10; 
<2 'Mean zenith angle and offset averaged over all contributing runs; 
<3) Acceptance corrected livetime of all H.E.S.S. observation runs pass- 
ing quality cuts; ( *'GC is detected with FermifLAT. 



such as position, distance (d), core radius (r c ), and central lumi- 
nosity density (po). To select the sample of GCs and observa- 
tional data that are suited for our study with H.E.S.S. we applied 
the following a priori cuts on the target and observation run list: 

1. GC IGalactic latitudel > 1?0 

2. Run passes standard quality selection criteria and is pointed 
<2°.0 offset from GC position 

3. GC has at least 20 available runs passing cut 2 

The first cut is to conservatively exclude GCs towards the di- 
rection of the Galactic plane to avoid chance coincidences with 
unrelated sources and to prevent potential contamination from 
a population of faint unresolved sources and/or diffuse emis- 
sion. Following the same approach as in the case of Terzan 5 
(see |Abramowski et al.|201 Ta) > we find that the probability of an 
individual GC coinciding by chance with an unrelated Galactic 
source based on its latitude is less than 10~ 3 for all GCs in this 
sample. The second cut excludes GCs that were only observed 
at very large offset angles from the camera center, because such 
observation runs would only provide relatively low effective ex- 
posure. With the third cut we make sure that there is a reason- 
ably long exposure on these potentially faint sources. In addi- 
tion we excluded Terzan 5 from this study because those results 
were recently covered in a separate paper (see Abramowski et al. 
201 la| ). Table[T]lists the 15 GCs passing these cuts. We note that 



one of the GCs, M 54, is actually not Galactic but most likely be 
longs to Sagittarius Dwarf Elliptical Galaxy (Ib ata et al.|19 94). 

Most of the GC positions were either serendipitously cov- 
ered in the H.E.S.S. Galactic plane survey (see Aharonian et al. 
2006b; G ast et al.|2012| l or lie in the same FoV as other sources 
observed with H.E.S.S. Therefore, the offsets of the GCs with 
respect to the camera center show a larger run-by-run scatter 
and in many cases larger mean values (see Tab. [T| compared 
to dedicated source observations. However, with a maximum al- 
lowed offset of 2°0 for individual runs (i.e. "cut 2") the dataset is 
still in a range where the analysis is very well calibrated. Also, 
the H.E.S.S. observations have been carried out during a time 



03 




longitude 

Fig. 1. Scheme showing the orientation of the four offset test 
regions with respect to the on-source region, aligned to the 
Galactic coordinate system. The offsets are CPA in both direc- 
tions and are defined such that region "test 4" is always closer to 
the Galactic plane than region "test 2". 



span of about six years where the gains of the photo multipli- 
ers and the optical reflectivity of the mirrors changed. To correct 
for these effects we calibrate the energy scale of the instrument 
using the response to single muons (see Aharonian et al. 2006a). 



2.2. Analysis of individual GCs 

We applied the standard quality selection cuts on the H.E.S.S 
observational runs to reject data affected by bad weather con- 
ditions or poor instrument performance (see | Aharonian et al. 



2006a). For the analysis we used the model technique where the 
air showers are described by a semi-analytical model. The ex- 
pected camera images are then compared to the observational 
data based on a maximum likelihood method (see de Naurois & 
Rolland 2009). The model analysis yields an improved sensitiv- 
ity, particularly at lower energies, and a more efficient Gamma- 
Hadron separation compared to the standard Hillas approach. We 
used standard cuts for the model analysis which include a 60 
photo electron cut on the image size. With this configuration the 
model analysis features a PSF with a 68% containment radius of 
CP.07 on average. We cross-checked all results with an indepen- 
dent calibration using a Hillas analysis framework, which em- 
ploys a machine-learning algorithm based on Boosted Decision 
Trees ( |Ohm et al.|2009l >. 

To search for VHE y-ray emission from each GC we per- 
formed two different kinds of analyses where the size of the 
emission region was assumed to be either a) point-like (CP. 1 inte- 
gration radius) or b) extended while using the nominal GC posi- 
tion as the center of the extraction region (see Tab. [TJ. For case 
b) we determined the expected intrinsic size (with the effects of 
the PSF removed) of the source related to each GC based on the 
measured extent of the VHE y-ray source HESS J1747-248 de- 
tected from the direction ofTerzan5 (Abramowski et al. 2011a). 
The morphology of HESS J1747-248 was fitted with an asym- 
metrical 2-dimensional Gaussian with intrinsic major and minor 
axes of CP.16 and CP.03, respectively. Assuming that the physical 
size of the VHE y-ray source is the same among all GCs, we 
scaled the angular size of the major axis of HESS J 1747-248 
with the relative distance of each individual GC (as given by 
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Table 2. Analysis results and model predictions 



GC 




WON® 


AW 2) 


l/a (3 > 


sig.< 4 > 


r (5) 






^UL / ^IC;IR,opt,CMB ^ 


name 


(TeV) 


(counts) 




(o-) 


n 


(phcm 


- 2 s"') 






a) point-like source 


analysis 


















NGC 104 


0.72 


72 


941 


18.2 


2.6 


- 


1.9x10" 


12 


2.6x10"' 


2.1x10' 


NGC6388 


0.28 


180 


2365 


14.9 


1.6 


— 


1.5x10" 


12 


8.0xl0" 2 


1.6x10° 


NGC 7078 


0.40 


119 


1988 


15.0 


-1.2 


- 


7.2x10" 


13 


1.9x10"' 


2.1x10' 


Terzan 6 


0.28 


202 


8194 


42.0 


0.5 


- 


2.1x10" 


12 


7.3x10"' 


1.0x10° 


Terzan 10 


0.23 


lb 


24D5 


36.0 


0.9 


— 


2.9x10" 


12 


4.3x10"' 


2.7x10"' 


NGC 6715 


0.19 


159 


2361 


15.2 


0.3 


- 


9.3x10" 


13 


3.1x10"' 


1.3xl0 2 


NGC 362 


0.59 


18 


533 


33.0 


0.4 


— 


2.4x10" 


12 


3.9x10° 


1.8xl0 2 


Pal 6 


0.23 


363 


10810 


31.4 


1.0 




1.2x10" 


12 


1.3x10' 


1.1x10' 


NGC 6256 


0.23 


64 


1869 


27.4 


U.J 




3.2x10" 


12 


1.8x10' 


2.9x10' 


Djorg2 


0.28 


56 


2387 


^0 4 


\J.KJ 




8.4x10" 


13 


1.0x10' 


1.0x10' 


NGC 6749 


0.19 


84 


2633 


70 3 


-0.6 




1.4x10" 


12 


2.5x10' 


4.1x10' 


NGC 6144 


0.23 


63 


2196 


JU.O 


— 1.0 




1.4x10" 


12 


5 Svlf) 2 
J.OA1U 


l.lxlO 3 


NGC 288 


0.16 


647 


24148 


JO.J 


U.O 




5.3x10" 


13 


7 7v1f) 2 

A. 1 A 1U 


q ?vin 3 

.J.^A 1U 


HP1 


0.23 


67 


2771 


^4 ^ 


— 1 

1 .u 




1.5x10" 


12 


5.2xl0 2 


1.7xl0 2 


Terzan 9 


0.33 


89 


2556 


^1 7 






4.5x10" 


12 


2.6xl0 4 


9.0xl0 2 


b ) extended source analysis 


















NGC 104 




293 


2016 


7.4 


1.2 


0.22 


2.3x10" 


12 


2.3x10"' 


1.9x10' 


NGC 6388 




253 


2818 


12.9 


2.2 


0.11 


1.7x10" 


12 


9.2xl0" 2 


1.8x10° 


NGC 7078 




161 


2386 


14.0 


-0.7 


0.11 


1.1x10" 


12 


2.8x10"' 


3.1x10' 


Terzan 6 




304 


9802 


34.2 


1.0 


0.12 


2.4x10" 


12 


8.1x10"' 


1.2x10° 


Terzan 10 




218 


A 1 Q A 

4i34 


19.0 


0.0 


0.18 


3.6x10" 


12 


5.4x10"' 


3.4x10"' 


NGC 6715 




159 


2361 


15.2 


0.3 


* 


9.3x10" 


13 


3.1x10"' 


1.3xl0 2 


NGC 362 




30 


708 


25.6 


0.4 


0.13 


2.5x10" 


12 


4.0x10° 


1.8xl0 2 


Pal 6 




1148 


17631 


16.6 


2.5 


0.18 


2.1x10" 


12 


2.4x10' 


1.9x10' 


NGC 6256 




131 


2524 


20.4 


0.6 


0.13 


3.9x10" 


12 


2.1x10' 


3.5x10' 


Djorg2 




137 


3753 


24.8 


-1.2 


0.16 


9.7x10" 


13 


1.2x10' 


1.2x10' 


NGC 6749 




168 


3544 


20.7 


-0.3 


0.14 


2.1x10" 


12 


3.6x10' 


5.9x10' 


NGC 6144 




120 


2913 


23.9 


-0.2 


0.13 


2.5x10" 


12 


6.7xl0 2 


1.9xl0 3 


NGC 288 




1030 


30767 


30.7 


0.8 


0.13 


6.1x10" 


13 


3.1xl0 2 


3.7xl0 3 


HP1 




67 


2771 


34.3 


-1.6 


* 


1.5x10" 


12 


5.2xl0 2 


1.7xl0 2 


Terzan 9 




206 


3909 


18.8 


-0.1 


0.16 
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"'Energy threshold of the analysis, defined as the location of the peak in the distribution of reconstructed photon energies; <2) Total number of 
on- and o ff-counts; (3) Rati o between off- and on-exposure when applying the reflected background technique; (4) Detection significance (pre-trial) 
following 



Li&Ma 



i jl983f ; (5) Extraction radius used for the an alysi s assuming extended emission; a star(*) denotes the cases where the calculated 
intrinsic extent is negligible compared to the PSF (see sect. 2.2 1; (6) Photon flux upper limits (99% c.L, following Feldman & Cousins 1998 1 
assuming a power-law spectrum with an index of -2.5; (7) Ratio between the upper limit and the expected IC flux using either the GC starlight 
(GC), or the sum of the Galactic infrared (IR), Galactic optical (opt) and the CMB background as target photon fields (see sect.|3]l. 



H10) and Terzan 5 (d T s = 5.9 kpc, |Ferraro et al.|2009| l. We note 
here that HESS J1747-248 is not only extended but also offset 
by 4!0+l!9 from the center of Terzan 5. Should such offsets be 
a general feature of VHE y-ray emission from GCs this would 
also affect the regions for analysis b). However, the reason for 
this offset has not been attributed to any physical parameter of 
Terzan 5 and particularly its direction for other GCs is therefore 
not easily predictable. Keeping the regions centered on the nom- 
inal GC positions but increasing the radius to also account for a 
potential offset could increase the fraction of the emission cov- 
ered by the extraction region. However, this would also heavily 
increase the number of background counts, thereby decreasing 
the sensitivity. Due to this fact, and since the offset seen from 
Terzan 5 is much less significant than the intrinsic extent of the 
signal, we decided to not include the offset in the determination 
of the expected source size. To define the final extraction region 



for analysis b) we assumed a 2-dimensional gaussian morphol- 
ogy with the width of the expected intrinsic size, folded with the 
H.E.S.S. PSE The radius of the extraction region was then the 
68% containment radius of the folded distribution. These radii 
are listed in Tab. E] Note that for the two GCs NGC 6715 and 
HP 1 the calculated intrinsic size was negligible (< 1% of the 
size of the PSF) and we therefore considered the point-like anal- 
ysis results also for case b). 

For the analyses we used the reflected background technique 
(see |Berge et al.|2007| ), which ensures that the regions used for 
signal (ON) and background (OFF) extraction have the same ac- 
ceptance in the FoV of the camera. In Tab. [2] we present the re- 
sults for both analyses and for each individual GC. From none of 
the GCs, neither with analysis a) nor b), significant excess emis- 
sion is seen above the estimated background, as the individual 
significances for a detection, even pre-trials, are well below the 
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Table 3. Results for on-source and offset regions 
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-0.31 


1.3 
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Mean and rms of the detection significance distributions for the on- 
source and offset test regions. Also, the total significances obtained with 
the stacking analyses are quoted. 



threshold of 5 cr. The mean (rms) values of the significance dis- 
tributions are 0.30 (1.0) cr and 0.34 (1.1) cr for analysis a) and b), 
respectively. Thus, the significance distributions are reasonably 
compatible with only background fluctuations within statistical 
uncertainties. We derived upper limits for the photon flux above 
the energy threshold (99% confidence level, following Feldman 
|& Cousins|1998| l assuming a power-law spectrum for the photon 
flux 



' P h' 



(— ) 
\TeV7 



(1) 



with an index of T = -2.5, which we chose to allow for a com- 
parison with the results obtained for Terzan5 (Tts = -2.5 + 
0.3 sta t ± 0.2 sys , |Abramowski et al.|201 la]> . Varying the assumed 
spectral index between -2.0 and -3.0 changes the upper limit 
by only 10% to 15%. 

As a cross check, we performed analyses with identical pa- 
rameters for four test positions in the vicinity of each GC. These 
regions are offset by ±(fA in Galactic coordinates (see Fig.[T]i. At 
this distance in no case a contamination from the supposed GC 
sources would be expected. The purpose of these test positions 
is to generate a sample of source-free regions with a similar dis- 
tribution in Galactic coordinates as the GC positions while still 
using a dataset comparable to the on-source regions. Also, be- 
cause the region test 4 is always closer to the Galactic plane than 
the other regions, particularly test 2, one can study the effects 
that the Galactic latitude may have on this particular population 
study. In Tab.|3]we show the parameters of the significance distri- 
butions for the on-source and the four test regions. As is evident 
from these results, we did not detect a significant signal from any 
of these regions. Even though there appears to be a slight ten- 
dency towards negative significances for the test regions, the ab- 
solute mean values are very small and nothing substantial can be 
concluded from that apparent shift. To estimate the background 
for the test regions we excluded the on-source area. Therefore, 
the slight shift towards negative significances for test 1 to test 4 
cannot be explained by a faint signal related to the GCs, which 
could have otherwise led to an overestimation of the background. 

2.3. Stacking analysis 

Because no individual GC studied in this sample shows signif- 
icant VHE y-ray emission, we performed a stacking analysis to 
search for a population of faint emitters. To combine the anal- 
ysis results and count-statistics from the observation runs of all 
GCs we applied the same method that is used to combine the re- 
sults for analyses of individual sources: The total number of ON 
(OFF) counts is the sum over all ON (OFF) counts of the runs 



in the stack. We calculated the ratio between the ON and OFF 
exposure (a) as the weighted mean of a of all observation runs 
with the weight being the number of OFF counts per run. The 
energy-dependent effective area of the stack is calculated as the 
average over all runs, weighted with the livetime. 

The total GC stack has an acceptance-corrected livetime of 
195 hours of good quality data and features an energy thresh- 
old of 0.23 Te V. We performed the stacking analysis for both the 
point-like and extended source analyses described in the previ- 
ous section. The results are shown in Tab. [2] The total detection 
significances obtained with these stacking analyses are 1.6 cr and 
2.4 cr, respectively, which are also very well below the threshold 
of 5 cr required for a firm detection. Therefore, we calculated up- 
per limits on the photon flux for the full stack, again assuming a 
power-law spectrum with an index of F — -2.5 (see Tab. [2]). 

For comparison we also performed the stacking analyses 
on each of the four offset test regions as defined in the previ- 
ous section. The stacked significances obtained for these regions 
are systematically lower than for the on-source positions (see 
Tab. [3}. However, the differences are not sufficiently significant 
to draw any firm conclusions at this point. 



3. Discussion 

Even though neither any individual GC nor the stack as a whole 
revealed a significant y-ray signal, the resulting upper limits on 
the photon flux can still be used to investigate models based 
on IC emission from leptons accelerated by msPSRs. All state- 
ments below rely on the assumption that the size of the potential 
emission region is covered by either of the two analysis meth- 
ods we employed for this study. Given the angular resolution of 
H.E.S.S., a point-like source would basically mean that the lep- 
tons are contained within a <4' radius around the GC center. For 
the analyses searching for extended emission we assumed that 
the physical size of the emission region is the same among all 
GCs, using the intrinsic extent of HESS J1747-248 as a tem- 
plate. This would be justified if HESS J1747-248 is indeed re- 
lated to the msPSR/IC scenario, if the leptons can propagate out 
of the half-mass region (for some smaller GCs even out of the 
tidal radius), and if the diffusion and energy-loss mechanisms as 
well as the age of the lepton population are comparable among 
all GCs. 

Now, assuming that HESS J1747-248 is indeed related to 
Terzan 5 and that the msPSR/IC scenario is the dominant process 
for VHE y-ray emission, the IC flux should scale as: 



F IC °C NmsPSR xN ph xd 



(2) 



with the number of msPSRs A^ ms psR, the total number of avail- 
able target photons in the volume occupied by relativistic leptons 
for the IC process N p h and the distance to the GC d (Eq.[2| is es- 
sentially the same as Eq. 8 in Venter & de Jager 200?). This 
simple scaling relation is of course only valid if the mechanisms 
for lepton production are similar among all GCs and, particu- 
larly, if the lepton energy spectra and the acceleration efficien- 
cies of the msPSR populations are comparable. In this case the 
proportionality constant in Eq. [2] can be fixed using the flux of 
HESS J 1747 -248, which then allows to calculate the expected 
VHE y-ray flux for the other GCs. In addition to the parameters 
given in Eq. [2] also the energy density of the ambient magnetic 
field in the emission region influences the IC flux due to the com- 
peting synchrotron radiation losses. However, this only changes 
our predicted IC fluxes if the magnetic field strengths (Bcc) for 
the 15 GCs are different from the one in Terzan 5 (Bts)- In this 
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case the expected IC flux should be scaled by an additional factor 
of b \sI b gc- 

The total number of msPSRs is difficult to access directly in 
a homogeneous way for all 15 GCs in this study. This is mainly 
due to the fact that the various radio surveys feature varying or 
unknown degrees of completeness, which becomes a very prob- 
lematic issue for such a large sample of GCs. As an alternative to 
the number count of detected radio msPSRs the total HE y-rays 
flux measured with Fermi/LAH could also be used as a direct 
indicator for the msPSR population. However, the GeV mea- 
surements suffer from large statistical uncertainties and are only 
available for two GCs in this sample. 

Here, we rely on a more indirect approach to estimate AVsPSR 
by using the stellar encounter rate F e = Pq 5 xr^ which is strongly 
correlated to the number of neutron star X-ray binaries, the pro- 
genitors of msPSRs (see, e.g. Pooley & Hut 2006 and references 
therein). Here po and r c denote the central luminosity density and 
the core radius of the GC, respectively. It has also been shown 
that T e is correlated with the HE y-ray flux ( jAbdo et al.|2010b 



Hui et al.|2010 i. However, there might be systematic uncertain 
ties in using T e as a proxy for A^ ms psR because the efficiency 
for the formation of msPSRs also depends on other parameters 
such as the metallicity of the GC (see e.g. Ivanova et al.|2 008). 
Assuming that F e is indeed proportional to A^ ms psR and that the 
proportionality constant is the same for all GCs, Eq.|2]then reads: 



(3) 



To calculate Fic we used the values for po, d and r c (con- 
verted to units of pc using d) as given by H10. For N v h we 
used four different target photon fields: 1) the optical light from 
the GC stars, 2) the Galactic infrared and 3) the Galact ic op- 
tical background photon fields (as used in GALPROHM Porter 



& Str ong 2005) at the location of the GC (as given by H10), 
and finally 4) the cosmic microwave background. The number 
of photons in case 1) is approximately proportional to the V- 
band luminosity of the globular cluster, which we chose for N v h 
using the values given in H10. For cases 2), 3) and 4) N v h is pro- 
portional to the energy density of the respective radiation fields. 
To compare the expected IC flux with the upper limits we first 
calculated for each GC the VHE y-ray flux of HESS J1747-248 
above the respective energy threshold, using the published re- 
sults (Abramowski et al. 2011a): flux normalization at 1 TeV: 
(5.2±l.l)xl0~ 1;t cm" 2 s" 1 TeV~\ and power-law spectral index: 
-2.5 + 0.3 sta t ± 0.2 sys . We then used this flux to determine the 
proportionality constant in Eq. [3] which in turn yields the ex- 
pected IC flux for the respective GC. Here we considered two 
extreme cases, namely that the leptons are either confined within 
the core region where the GC starlight is the dominant target 
photon field, or that the particles have propagated far away from 
the core where the Galactic background fields and the cosmic 
microwave background become the most important components. 
In the latter case we used the sum of the energy densities of the 
three photon fields as a proxy for N p h. The ratio between the flux 
upper limits and the expected IC flux are given for both cases in 
the last two columns of Tab. [2] 

In the simple scaling approach (Eq.[3]l to estimate the IC flux 
we assumed that the lepton spectral shapes as well as the ac- 
celeration efficiencies of the msPSR populations are comparable 
among all GCs, which might not be the case. Potential deviations 
could arise from different energy-loss and diffusion timescales, 
differences in particle (re-)acceleration mechanisms or the spec- 
tral shapes of the injection spectra. Those and other issues would 
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need to be addressed by more detailed models for each individ- 
ual GC. Apart from the systematics mentioned above, the param- 
eters that are directly used in Eq.[3]are also subject to uncertain- 
ties, which are discussed below. The most impact by far for po- 
tential uncertainties comes from the distance d because Fie itself 
depends inverse-quadratically on d; also the values of po as well 
as r c are less strongly dependent on d. Unfortunately, the uncer- 
tainties of d are perhaps the most difficult to quantify because 
there are several different methods for distance determinations 
of GCs, each suffering from other systematic issues. Therefore, 
we decided here to use the distances quoted by H10 for all of the 
GCs studied in this work, to only rely on one method. Assuming 
an error margin for d of a factor of two (based on typical vari- 
ations between different distance estimation methods) yields an 
uncertainty of a about a factor of five in the expected IC flux. 
Additionally, the confidence intervals of the spectral parameters 
of HESS J1747-248 ( |Abramowski et al.||2011a| ) introduce fur- 
ther uncertainties in the expected IC flux of about 20 %, yielding 
a total error margin of a factor of ~5. To estimate the uncertain- 
ties of the predicted model flux from the full stacks, we assumed 
a factor of 5 uncertainty for each individual GC and propagated 
the errors. We show the resulting uncertainty band for the stack- 
ing analyses in the last two columns of Tab. [2] 

With the optical light of the GC member stars as the tar- 
get photon field, the flux upper limits for point-like emission 
from the two GCs NGC6388 and NGC7078 lie about one or- 
der of magnitude below the expected IC flux (see Tab. [2]). This 
is mainly due to the fact that Terzan 5 is relatively faint in the 
V-band compared to many other GCs, which yields large ex- 
pected IC fluxes for GCs with a similar or even larger stellar 
encounter rate F e . Among all GCs studied here, given the un- 
certainties related to the predictions, only in these two cases the 
simple scaling model is challenged by our observational results. 
However, the fact that HESS J1747 -248 is extended and offset 
beyond the half-mass-radius of Terzan 5 could point towards fast 
diffusion and long energy-loss timescales of the lepton popula- 
tion in Terzan 5. If this is indeed the case, then other target pho- 
ton fields would be dominant over the GC starlight. Therefore, 
the results of the point-like analyses are better suited to be com- 
pared to the predictions with the GC starlight as target photon 
fields, whereas the results from the extended analyses are more 
relevant for the predictions based on Galactic background pho- 
tons. We note here that for some GCs the Galactic background 
could still make a contribution to the IC flux even for regions 
close to the GC cores, an issue that is not taken into account in 
our scaling model. Looking at the results for the infrared and op- 
tical background photon fields, the upper limits do not allow for 
similarly strong constraints as in some cases when considering 
the GC starlight. Given the uncertainties related to the parame- 
ters in Eq. [3] the msPSR/IC scaling model cannot be ruled out 
for any of the GCs assuming Galactic background photons as IC 
target. 

To make a statement about the GC sample as a whole, we 
derived the expected IC flux from the total stack by calculat- 
ing the weighted mean flux from all GCs. As weights we used 
the acceptance-corrected livetimes of the individual clusters. The 
upper limits for the stack are about two orders of magnitude 
below the expectation for the GC starlight and at a level of 
~50% for the external target photon fields (see bottom section of 
Tab. [2|. In the first case, even considering the uncertainties, this 
seems to be quite a strong constraint on the simple msPSR/IC 
scaling model. Therefore, either the IC flux from GCs does not 
scale as expected, or the VHE y-ray source near Terzan 5 is not 
related to msPSRs. However, for the second case, i.e. assum- 
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ing background photon fields as the dominant IC target, the up- 
per limits for the stack are compatible with expected fluxes, and 
therefore do not allow to put similarly strong constraints on the 
model. 

There are many possibilities why Terzan5 may stand out 
among the general GC population. In the msPSR/IC model the 
confinement of leptons in the emission region might be stronger 
in Terzan 5 compared to other GCs, which would lead to an ac- 
cumulation of particles from earlier epochs that would there- 
fore enhance the observed IC flux per msPSR. However, be- 
cause HESS J1747-248 is extended and offset from the core, 
the trapping mechanism is required to work on similarly large 
spatial scales. A high interstellar magnetic field strength in the 
vicinity of Terzan 5 could provide such a trapping mechanism, 
but would also increase the synchrotron losses and therefore de- 
crease the expected IC flux. These competing effects would need 
to be modeled in detail for each individual GC, and the results 
presented in this paper might yield interesting constraints. 

As already mentioned in sect. [T] there are also models that 
are not related to the msPSR/IC scenario which suggest that 
the VHE y-ray emission is due to hadronic processes in a short 
gamma-ray burst remnant (Domainko]|201 lb. Due to the ex- 



tremely low rates of these catastrophic events of only -10 yr 



(Nakar 2007| ) it could be easily explained why only one GC 
stands out among the whole population. 

4. Conclusions 

Judging from recent HE and VHE y-ray observations, GCs 
could feature interesting high-energy processes, and, in partic- 
ular, could constitute a new class of sources in the VHE y-ray 
regime, if the recently detected source HESS J 1747-248 is in- 
deed related to the GC Terzan 5. 

To test the validity of the msPSR/IC scenario for VHE y- 
ray emission, we analyzed the H.E.S.S. data of 15 Galactic GCs. 
According to our observations, neither any individual GC nor 
the stacking analyses show significant excess emission above the 
estimated background. Based on a simple model for IC emis- 
sion and the measured VHE y-ray flux of HESS J 1747 -248 we 
calculated predictions for the fluxes for each individual GC and 
compared them with the derived upper limits. 

The non-detections of the two GCs NGC6388 and 
NGC7078 challenge the simple IC scaling model when as- 
suming the GC starlight as the dominant target photon field. 
However, when searching for VHE y-ray emission from a larger 
scale, comparable to what is seen from Terzan 5, and consid- 
ering Galactic optical, infrared and CMB photons as the dom- 
inant IC target, no clear constraints can be put on our simple 
msPSR/IC scaling model. The upper limits derived from the 
stacking analyses are about two orders of magnitude below the 
expected cumulative IC flux from the whole GC sample when 
assuming the GC starlight as the dominant target photon field. 
Taking the uncertainties of those parameters into account that 
directly influence our flux predictions, these results still exclude 
the simple msPSR/IC scaling model for this case. But, as men- 
tioned in sect. [3] there are several caveats to the simple scaling 
model applied here that would have to be considered in more 
elaborate and custom-tailored models for each individual GC. 
Also, there is still the possibility that the VHE y-ray signal from 
Terzan 5 might not be related to the msPSRs at all and arises 
from hadronic or other processes. 

Further studies of GCs are still ongoing and future deep ob- 
servations with Cherenkov telescopes and X-ray observatories of 
the most promising candidates will increase the strength of the 



derived upper limits, or might even facilitate another detection 
next to Terzan 5. 
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